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Abstract 

We calculate the exact dynamical magnetic structure factor S(Q, E) in the 
ground state of a one-dimensional S = 1/2 antiferromagnet with gapless 
free 5 = 1/2 spinon excitations, the Haldane-Shastry model with inverse- 
square exchange, which is in the same low-energy universality class as Bethe's 
nearest-neighbor exchange model. Only two-spinon excited states contribute, 

and S(Q, E) is found to be a very simple integral over these states. 
75.10.Jm, 71.10,+x 
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The integrable "Haldane-Shastry" model |1]|| (HSM) is a variant of the S = 1/2 Heisen- 
berg chain, with exchange inversely proportional to the square of the distance between spins, 
in contrast to the nearest-neighbor exchange of the Bethe Ansatz model || (BAM). The 
low-energy properties of the antiferromagnetic HSM are in the same universality class as 
those of the BAM: the elementary excitations are S = 1/2 objects || ("spinons") obeying 
semion (half-fractional) statistics intermediate between bosons and fermions f|, and the 
low-energy fixed point is described by the k = 1 SU(2) Wess-Zumino-Witten model, a c = 1 
conformal field theory |l| . A generic model in this universality class has marginally irrel- 
evant interactions that renormalize to zero at the infra-red fixed-point; this is seen in the 
BAM. The special feature of the HSM is that such interactions are absent at all energy 
scales ||: it is in a very real sense the model in which the spinons form an "ideal S = 1/2 
semion gas" , and may be regarded as the fundamental model for gapless half-integral spin 
antiferromagnetic chains. 

While the correlation functions of the BAM have not yet been obtained, the static 
(equal time) two-point antiferromagnetic ground-state spin-correlations ofthe HSM are al- 
ready known JT|fJ] . In this Letter, we extend this result to the full dynamical spin correlation 
function: 

(0\S a m (t)S b n (t')\0) = \5 a \-lT- n C{m-n,t-t') (1) 

where our expression for C(x,t) is remarkably simple, and directly related to the spinon 
spectrum. We will first present the result, then describe the derivation. 
The HSM Hamiltonian is 

H = JY,(d(m-n))- 2 S m -S n , (2) 

where d(n) = (N/ir) sin(7m/iV) —>■ n as the number of sites N — > oo. The HSM spinon 
dispersion relation is 

e(q) = (v/n) ((vr/2) 2 - q 2 ) , v = vr J/2, (3) 
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which is restricted to half the Brillouin zone (\q\ < vr/2); v is the low-energy spinon velocity. 
For iV — > oo, a state with N sp spinons has energy Y^,i e (%) and crystal momentum exp(iK) = 
(— 1) M exp(z J2i Qi), where M = (N — N sp )/2 must be integral. 
Our concise expression for C(x,t) is 

J dM exp (±ivrTr[xMM - \vt{l - (MM ) 2 )]) . (4) 

The integral is over the manifold of 4 x 4 traceless Hermitian unitary matrices M = M _1 
= (&l • (Tl){Q,ii ■ (Tr), where {cr£} and {&%} are two independent sets of 4 x 4 Hermitian 
generators that individually obey the same algebra as the Pauli matrices, and which commute 
with each other; Cl i and fi^ are real unit 3- vectors. This manifold is isomorphic to the 
product of two spherical surfaces S\ x S 2 R . The time-reversal operator for Pauli matrices 
is r = icr 2 , so (cr a )* = T<r a T, and r 2 = — 1; M obeys the reality condition M* = (t^t^) 
M(tlt#) where (tlTr) 2 = 1. Mo = (z • «?l)(z • <t#) is a point on the manifold. Tr[MMo] 
= Ailing and Tr[l - (MM ) 2 ] = 8((fil) 2 + (fifj 2 - 2(fi^) 2 ). 

The invariant measure for the integral is the product of rotationally-invariant measures 
on S 2 x Sft, and the normalization is fixed so C(0, 0) = 1. Writing fi z L ^ = Ai( 2 ), this gives 
C(x, t) as 



1 1 1 tZAi f d\ 2 e iQx ~ Et (5) 

1 J — 1 
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where Q = 7rAiA 2 and E = (7rw/2)(A 2 + A 2 — 2A 2 A 2 ). Alternatively, if 7rAiA 2 = qi + (fe, and 
(Qi ~ Q2) 2 = tt 2 (1 - A 2 )(l - Af>), so Q = gi + g 2 and E = e{q x ) +e(q 2 ), C{x,t) is also given by 



(2tt) 2 J-tt/2 y-Tr/2 [e{qi)e{q 2 )\ 1/2 

Note that the 5* = 1 state S^|0) is completely expressible in terms of eigenstates of the HSM 
with only two parallel-spin spinons carrying momenta q x and q 2 . Thus (§) explicitly expresses 
the correlation function in terms of the physical excitations; this formula is the principal 
result reported here. 

(0|5^(t)$£(f )|0) can be expressed in terms of the dynamical structure factor S(Q, E) as 



' dQ dES(Q,E)e i ^ m - n ^ E ^\ (7) 



(8) 



2n JO JO 

where S(Q, E) is explicitly given by 

~e{E-E lv {Q)) 



\9{E 2 {Q)-E) J] 



[E-E lv {Q)] l/ 



where E X _{Q) = {v/u)Q{u - Q), E 1+ {Q) = (v/ir)(Q - n)(2n - Q), and E 2 (Q) = 
(v/27t)Q(2tt — Q). Integration over E gives the previously-known expression for the static 
structure factor: S(Q) = |log(|l — Q/n^ 1 ) . 

Our calculation is based on the identification of C(x,t) as proportional to the bosonic 
single-particle ground-state correlation function (0|\l/t(a/, £')v|/(0, 0)|0) (with x' oc x, t' oc t) 
for the Calogero- Sutherland model (CSM) of a spinless ID Galilean-invariant gas with 
interactions (U 2 / Am)((3(f3 — 2)/(x i — Xj) 2 ), at coupling constant @ — 4, the coupling at which 
it is related to the symplectic random matrix ensemble |J. 

Though different in its technical details, a recent calculation by Simons, Lee and Alt- 
shuler (SLA) of (0\p(x', f)p(0, 0)|0) for the (3 = 4 CSM suggested our calculation. Our 
result shows that the hole excitation in this model fractionalizes into two semions (as in the 
essentially similar bosonic Laughlin fractional quantum Hall effect (FQHE) state at Landau 
level filling v = 1/2). Our interpretation of the SLA expression || is that it shows that 
the particle excitation remains elementary (the FQHE analog is a particle added outside 
the Laughlin droplet); inspection of their formula shows it to be the convolution of our re- 
sult for the fractionalized hole contribution with a factor representing an additional particle 
excitation outside the ground state condensate. 

We first note that (for finite N) the HSM exhibits not just SU(2) symmetry, with an 
associated sl 2 Lie algebra generated by Jo = J2m S m) but a remarkable additional quantum 
symmetry, the Yangian Y{sl2) symmetry algebra with the additional generator [|Kj 



h \— v (itim — n)\ _ „ . . 

J i = o E cot Ar S m x S n (9) 
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where [H,J^] = J^JO) = 0, m = 0,1. Here h is the "quantum parameter" of Yfafa); this 
is a scale parameter, conventionally rescaled to 1. However, in the limit N —>■ oo, hN 



must be held constant, so h — > 0. The symmetry algebra then is the "classical" infinite- 
dimensional Lie algebra SI2+, which is the algebra of non- negative modes of a Kac- Moody 
algebra: [J£, J b n ] = ie abc J^ +n , m,n > 0, with J£|0) = 0, m > 0. 

HSM eigenstates are organized into Yangian "supermultiplets" (of many degenerate 

577(2) multiplets) each containing a single Yangian highest weight state (YHWS) which is 
the state of highest Jq in the multiplet. The set of YHWS span a subspace of states 



that can be physically characterized as the "fully-polarized spinon gas" (FPSG) states [10 
where all spinons have maximally polarized spins. The wavefunctions for FPSG states can 
be conveniently expressed in terms of M complex coordinates Zj (with (Zi) N = 1) of lattice 
sites (on the unit circle) with reversed spins: 

*{{Zi}) = U(Zi - z 3 ) 2 Y[z tJ (io) 

i<j i 

where $({Zj}) is a symmetric polynomial of degree N sp in each variable; the constraint 
M = (N — N sp )/2 ensures that the Fourier transform of the wavefunction has compact 
support within a single Brillouin zone. Diagonalization of the HSM Hamiltonian within the 
FPSG subspace gives the YHWS. The other eigenstates are obtained from the YHWS by 
repeated action of the Yangian generators, and the multiplets are classified as irreducible 



representations of Y^sfa) [^0| , pj| . The structure of these multiplets is reminiscent of that 
of the ideal S — 1/2 Fermi gas, where there is an independent spin degree of freedom for 
each orbital; Y(s/ 2 ) multiplets are also isomorphic to direct products of independent 
representations, but in contrast to the fermion case, the "orbital" or "Fock space" structure 
changes with particle number. These degeneracies make it appropriate to identify the 
excitations of the HSM as forming an "ideal spinon gas" [|B , 10f . 



The operator £"+ acting on a FPSG state removes a down-spin coordinate, and preserves 
the FPSG character of the state. The wavefunction describing S'+|0) is 

U.{Zi-*n) 2 U(Zi- z i) 2 U z i- ( n ) 

i i<j i 

The fact that the initial polynomial factor is degree-2 in each coordinate Zj immediately 
shows that this state is composed only of states with two spinon excitations. 



The polynomial character of the FPSG wavefunctions implies that when matrix elements 
of the Hamiltonian between such states are calculated, sums over discrete lattice sites on the 
unit circle can be replaced JTJQ by integrals of a continuous coordinate on the unit circle (the 
compact support of the Fourier transforms of these polynomials within the first Brillouin 
zone in reciprocal space means that all "Umklapp" corrections vanish). Thus calculations 
for the HSM involving only FPSG states will give identical results to analogous calculations 
for a continuous model, the A = 2 CSM on a circle. Then C(x,t) oc (0\^(x', t')^(0, 0')|0) 
where x' and t' are space and time coordinates suitably scaled in terms of CSM units. 

Because of the special form of the HSM and CSM ground states, C(x,t) can also be 
written 

<0| U(Z t - z'fe-^ H - E ^' H(Z t - ^) 2 |0). (12) 

i i 

We note that the method of our calculation can easily be extended to obtain 
(0| Ylj W(xj, t) Ylj ty(x'j, t')\0) where n particles are removed at various places at one time t' 
and replaced at different places at a different time t. (The final integration over a manifold 
of 4n-dimensional matrices will be more complicated, however.) 

In the thermodynamic limit, the correlations of the CSM for particles confined in a 
weak harmonic potential that is rescaled to maintain a constant particle density in some 
region become locally equivalent to those of the same density CSM model on a circle. We 
use this Gaussian formulation, since in this case ||12|| , the dynamics of the CSM particles 
become identical to the dynamics of the doubly-degenerate eigenvalues of a Hermitian matrix 
with matrix elements that are harmonic oscillators constrained so the matrix has Kramers 



degeneracy. This is the Gaussian Hermitian matrix-model [12| with symplectic symmetry. 
The action for the matrix-model is 

S = I J dt Ti[(Q) 2 -lu 2 (Q) 2 ] (13) 

where Q is a 2N x 2N Hermitian matrix 

Q = S® l + ^A a ® a a . (14) 



Here S is a real symmetric N x N matrix and A° are three independent imaginary antisym- 
metric matrices, making 2N 2 — N independent coordinates in total. Then Q* = TQT -1 , 
where T = 1 eg) t, and T 2 = —1, which ensures Kramers degeneracy. The limit N — > oo 
with uj oc 1/N is taken at the end of the calculation. 

The SLA calculation || was originally a quenched random matrix ensemble average, but 
can be reformulated |13] in terms of a dynamical Gaussian Hermitian matrix-model, when 



it is equivalent to a calculation of 



(0\T t J] Tr[(Q(ai) -atx± iO H 

a=± 



,-1 



which is a more complex calculation than the one we describe here. 

In fljjP, if the {Zi} are interpreted as the doubly-degenerate eigenvalues of a unitary 
matrix U with Kramers degeneracy belonging to Dyson's circular ensemble with symplectic 
symmetry ||, Hi(Zi — z) 2 = det[U — z] . Invoking the large- N local equivalence of the 
eigenvalue distributions of the circular unitary and Gaussian Hermitian ensembles |J, we 
identify the Gaussian matrix-model quantity we need as 

(0\T t Y[ det[Q(at) - ax]|0) (15) 

a=± 

where the space and time coordinates have been rescaled in terms of matrix model units, 
and will be held fixed as iV — >• oo. 

The fundamental matrix-model correlation function is 

(0|T,Tr[AQ(0]Tr[BQ(O]|0) = G(A, B)g(t - t') (16) 

where g(t) = (h/2u) exp(iu\t\), A and B are c-number matrices, and G(A, B) is 

I E (^f B t + A u n ''/ - A if B fj a ) ■ ( 17 ) 

ijaa' 

The determinant is given by a Gaussian Grassmann-number integral: 

det[Q(at)] = / JJ d^^d^a) expT>[A(a)Q(a*)] 



where A°f (a) = ijj ia (a)ip* a ,(a). 

We now evaluate the harmonic oscillator correlation function, using the standard har- 
monic oscillator result (0|T t exp Ai exp A 2 \0) = exp[(l/2) J2ij{^\TtAiAj\0)} where Ai are lin- 
ear functions of the coordinates, and obtain 

I] / diPl(a)diPUa)expS[{^l(a)^Ua)}]; (18) 

iaa 

here S = (xSi + (Tl/8uj)S 2 ) where Si = Y,a a Paa and S 2 = Y,ai3(9ai3(t)) 2 G(A(a) , A(/3)) , with 

G(A(a), A(/3)) = - ( Pa(3Pf3a + \ct>* a ^ Pa ) , (19) 

and g a p(t) = 5 a/3 + (1 - 5 a p) exp(iu\t\). Here p a p = E iCT V4(a)^(/3), 0^/3 = 
Ew T C(jl il) i(J (a)il) ia ,((3) and 0*^ = £w T^^ia)^,^). 

We now carry out a bosonic Hubbard-Stratonovitch transformation to decouple S 2 ; ten 
independent real integration variables are required, which can be organized into a 2 x 2 
Hermitian matrix V a p coupling to p a/ 3 and a 2 x 2 complex symmetric matrix A Q| g coupling 
to (p al 3, with A*^ coupling to 0*^. These can in turn be organized into a 4 x 4 matrix 
M(V, A, A*) where M al , pi = V af3 , M al)/32 = A a/3 , M q2i/ ji = A*^ and M a2 ^ 2 = Vp a . Provided 
Vpa = (V a p)* and A*^ = (A^)*, this matrix is Hermitian, and obeys the reality condition 

M* = (1 <g> cr 1 )M(l <g> cr 1 ). (20) 

Such a matrix has the ten-parameter form 

M = C°l + ^2 G ah cr a L cr b R . (21) 

ab 

where C° is a real number and C is a real 3x3 matrix. The explicit realization of the 
generators is cf L = (cr 2 <g> cr 2 , — cr 2 <g> cr 1 , 1 <g> cr 3 ) and cf R = (cr 1 (g) cr 3 , cr 2 <g> 1, cr 3 <g> cr 3 ). 
The unitary transformations generated by {cr a L , cr a R } correspond to a G = SO(3)l® SO(3)r 
continuous symmetry group, under the action of which C — > O^CO^ 1 . 

We now take the large- N limit with 00 oc 1/N, and x',t' ~ 0(1). The Grassmann 
integrals decouple into a product of N identical integrals leading to a factor 
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(det[M + iN^SM]) 



N 



(22) 



where (dropping scale factors in x' and t') 



5M= (x' - t'[M,M ])M ( 



(23) 



Mo = (er 3 <g> 1). The correction <5M shifts the matrix in the determinant resulting from the 
Grassmann integrals away from the Hermitian value M of the Hubbard-Stratonovich fields, 
but this does not affect the validity or convergence of the Grassmann integral. Using the 
identity det[M] = exp(Tr[lnM]), the value of the integrand for large N, with the Gaussian 
Hubbard-Stratonovich factor included, becomes (det[M] exp(— Tr[M 2 ]/2)) 7V exp(i<p) where 
if = Tr[M _1 (x' — £'[M, Mo])Mo] • The factor which is exponentiated to a power N is 
maximized when M 2 = 1, or M _1 = M, so at the extremal points, 



There are five classes of solutions of the condition M 2 = 1 corresponding to the five 
possible signatures s = ±4, ±2, and of M, all with Tr[M] = s and C° = s/4. For s = ±4, 
there are two discrete solutions with C = 0, invariant under the full group G. For s = ±2, 
C = C O, where O is orthogonal, which is invariant under the SO (3) subgroup of G where 
Ol = OOrO~ 1 , giving two manifolds of dimension 3. Finally, for s = 0, C ab = fiJ^R where 
&l(r) are un it 3-vectors; this is invariant under the 0(2) L ® 0(2) R subgroup of G where 
Ol(r) is a rotation around flum, and forms a manifold of dimension 4. 

When the same extremal value occurs on several different manifolds, only manifolds with 
the largest dimension contribute to the integral as iV — > oo. The measure for the integral 
over the dominant s = extremal manifold is just the invariant measure on Si x S R . To 
complete the derivation of (|J), the space and time variables x' and t' must be rescaled into 
HSM units. This is accomplished by noting that the two-spinon states span a momentum 
range of 2ir and an energy range of nv/2. 

Finally, we comment on the relation between our HSM result and the BAM corre- 
lations. An Ansatz similar to the exact HSM result (||) for S(Q, E) has previously been 



<p = Tr[x'MM - £'(1 - (MM ) 2 )]. 



(24) 
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proposed |14| as an approximation to the BAM structure factor. The key simplification in 
the HSM is that only the two-spinon excited states contribute to S(Q, E) at T — 0. In the 



BAM, the numerical studies [II] show that while these states dominate the spectral weight, 
there is a finite contribution from the N sp > 2 states which, while small, presumably remains 
finite for any large but finite N sp . 

A variant of the HSM with N = oo and d(n) = k _1 sinh(/tn), k real, has been identified 
1 15|] as a family of integrable models that interpolate between the N = oo HSM (k = 0) 



and the BAM (k = oo). These models also have a Y(sl2) quantum symmetry [T(| generated 
by the analog of (H) with cot(7r(m — ri)/N) replaced by icoth(fi;(m — n); this realization of 
Y(sl2) goes over smoothly into its more-familiar realization in the N = oo BAM limit. 

The Yangian quantum parameter h must be rescaled to zero to keep h/n constant in the 
HSM limit k — > with iV = oo. The reverse process by which the "classical" 3/2+ symmetry 
at k = becomes the Y(sl2) "quantum" symmetry when k > is a "quantum deformation" . 
We speculate that the extension of the k = result (|6]) to the k > model |15| , and thence 
to the k = 00 BAM limit, may be achievable as such a "quantum deformation". In this 
scenario, there is an expression for C(m,t;K,) as an expansion in multi-spinon terms with 
N sp = 2, 4, 6, . . . , 00; when k = m /N the N sp > 2 terms would vanish, giving the finite- N 
HSM correlations. 

This work was conceived and carried out at the Aspen Center for Physics during its 1993 
Summer Program. MRZ acknowledges partial support from the Sonderforschungsbereich 
341, Koln-Aachen-Julich. 
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